Nucleon-Nucleon Interactions from Dispersion Relations: Coupled Partial Waves 
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We consider nucleon-nucleon interactions from chiral effective field theory applying the N/D method. The 
case of coupled partial waves is now treated, extending Ref. [1], where the uncoupled case was studied. As 
a result, three N/D elastic-like equations have to be solved for every set of three independent coupled partial 
waves. As in the previous reference the input for this method is the discontinuity along the left-hand cut of 
the nucleon-nucleon partial wave amplitudes. It can be calculated perturbatively in chiral perturbation theory 
because it involves only irreducible two-nucleon intermediate states. We apply here our method to the leading- 
order result consisting of one-pion exchange as the source for the discontinuity along the left-hand cut. The 
linear integral equations for the N/D method must be solved in the presence of I — 1 constraints, with I the 
orbital angular momentum, in order to satisfy the proper threshold behavior for £ 3> 2. We dedicate special 
attention to satisfy the requirements of unitarity in coupled channels. We also focus on the specific issue of 
the deuteron pole position in the 3 S \- 3 Di scattering. Our final amplitudes are based on dispersion relations and 
chiral effective field theory, involving only convergent integrals. They are amenable to a systematic improvement 
order by order in the chiral expansion. 
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I. INTRODUCTION 

Recently we employed, in Ref. [1], the N/D method [2] 
to study nucleon-nucleon (AW) uncoupled partial-waves in 
connection with chiral perturbation theory (ChPT) [3, 4]. In 
this approach the chiral counting is applied to the imaginary 
part of the AW partial waves along the left-hand cut (LHC), 
owing to (multi-)pion exchanges. This can be done because 
Cutkosky's rules require putting pion lines on-shell to cal- 
culate the discontinuity across the LHC, giving rise to irre- 
ducible nucleon diagrams. For more details see Ref. [5]. At 
this point we avoid calculating perturbatively contributions 
that involve A^-nucleon reducible graphs. 1 This method pro- 
vides AW partial waves by solving a linear integral equation 
that, by construction, involves convergent integrals and sub- 
traction constants that can be calculated in terms of physical 
quantities. As a result, no need for any type of cutoff arises in 
our novel approach. 

The idea of applying ChPT to evaluate irreducible N- 
nucleon contributions was originally put forward in Refs. [6]. 
There it was applied to calculate the effective multi-nucleon 
potential, which is later implemented in a Lippmann- 
Schwinger (LS) equation (or Schrodinger equation) in order 
to derive the full S -matrix. However, owing to the singular 
nature of the chiral potentials resulting from their calculation 
in ChPT the solution of the LS equation requires of some kind 
of regularization, typically a three-momentum cutoff A [7-9]. 
Several works [10-14] have shown that the chiral countert- 
erms that appear in the ChPT potential following the standard 
ChPT counting [6] are not enough to reabsorb the cutoff de- 



pendence that stems from the solution of the LS equation. Sta- 
ble results with the AW potential determined from one-pion 
exchange (OPE) are obtained in Refs. [10, 15] for A — > oo, by 
promoting counterterms from higher orders to lower ones. 2 
This implies a violation of the standard ChPT counting and of 
the low-energy theorems relating the parameters in the effec- 
tive range expansion [16]. One counterterm is needed for each 
partial wave with an attractive OPE tensor force [10], so that 
the AW scattering amplitude from the OPE potential would re- 
quire an infinity number of them, so that it is nonrenormaliz- 
able. A result compatible with this conclusion is also obtained 
in Ref. [17], where it is found that the AW scattering ampli- 
tude from the OPE potential is nonrenormalizable unless the 
tensor force part vanishes for A — » oo [17]. One should be 
aware that when A — > oo a more involved counting emerges 
[17-19]. The extension of these ideas to higher orders in the 
chiral potential is not straightforward and, up to now, cannot 
avoid cutoff dependence [15, 20]. On the other hand, the ap- 
plication of Weinberg's scheme has given rise to a great phe- 
nomenological success in the reproduction of AW phase shifts 
if the cutoff is fine-tuned in a region around 600 MeV, not be- 
yond the breakdown scale of the effective field theory (EFT) 
[8, 9]. Of course, the cutoff dependence is not removed then. 

We present the generalization of Ref. [1] to the case of 
coupled channels in Sec. II, where the corresponding three 
linear integral equations needed for each set of coupled par- 
tial waves are derived. We apply this method to leading order 
(LO), which implies taking OPE as the source for the discon- 
tinuity along the LHC. References [21-24] applied the N/D 
method to study AW scattering quantitatively. Reference [21] 
was restricted to the S waves and took only OPE as input 
along the LHC. References [22, 23] included other heavier 
mesons as the source for the discontinuity along the LHC, in 
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1 The latter require an extended version of the standard chiral counting as 
derived in Ref. [5]. 



In Ref. [10] the cutoff range was taken A < 4 GeV. 



2 



line with the meson theory of nuclear forces, so popular those 
days, while Ref. [24] modeled the LHC discontinuity by OPE 
and one or two ad-hoc poles. We stress that we present here a 
novel way to introduce the N/D method in harmony with the 
modern perspective of EFT. In this way, we show that one can 
calculate systematically within ChPT, according to the stan- 
dard chiral counting, the discontinuity along the LHC that is 
the basic input for the N/D method. This allows one to im- 
prove the results order by order, which was not the case by 
applying previous schemes [24-27]. This is a point of fore- 
most importance. 

In addition, the threshold behavior of partial waves with or- 
bital angular momentum 6 > 2 is satisfied within our approach 
by including zeroes at oo in the NN partial waves, which is al- 
ways allowed within the N/D method (these are the so called 
Castillejo-Dalitz-Dyson (CDD) poles [28]). However, in the 
previous works [22, 23] the correct threshold behavior was 
achieved in an ad-hoc way by including a fictitious pole below 
threshold, with the subsequent dependence of the results on its 
location, which was fitted to data. Furthermore, in our results 
we always respect coupled-channel unitarity, which was not 
the case in Refs. [21-24]. 

The results obtained with this formalism are first consid- 
ered for 3 S i- 3 £>i coupled waves in Sec. Ill A, where the spe- 
cific issue of the deuteron pole is also discussed. Higher par- 
tial waves are considered in Sec. IIIC. Our conclusions are 
collected in Sec. IV. Finally, we show in the Appendix the 
cancellation of a potential divergence in a function involved 
in our equations. This cancellation occurs thanks to the con- 
straints already imposed to satisfy the right threshold behavior 
for partial waves with £ > 2. 



II. COUPLED PARTIAL WAVES 

For spin triplet AW partial waves with total angular mo- 
mentum / one has the mixing of the orbital angular momenta 
i = J - 1 and t' — J + 1 (except for the 3 Po partial wave). 
Each set of coupled partial waves is determined by the quan- 
tum numbers S , J, i, and €' , where S is the total spin. In the 
following, to simplify the notation we omit them and indicate 
the different partial waves by fy, with i = 1 corresponding to 
t = J — 1 and i — 2 to (' — J + 1 , a convention that we adopt 
henceforth. As a result, a two-coupled-channel T-matrix re- 
sults. In our normalization, the resulting S -matrix reads 



S(A) = I + i2p(A)T(A) = 



cos2ee 2i<5 > ;sin2ee' x * 1+<52) 
/sin2ee /(,5l+t52) cos2ee 2,<52 



(1) 



such that <5i corresponds to the phase shifts for the channel 
with { = J - 1 and 6 2 to that with V = J + 1 . The argument A 
refers to the center-of-mass (CM) three-momentum squared. 
We have indicated the phase space by p(A), which reads, in 
our non-relativistic approximation, p(A) — m^A/An, where 
m is the nucleon mass. 

Along the right-hand cut (RHC), which corresponds to the 
physical region with A > 0, the unitarity of the S -matrix, 



S S 1 = S T S - I, can be written in terms of the (symmetric) T- 
matrix as Imr -1 (A) = -p(A) I. In the following, the imaginary 
parts above threshold of the inverse of the T-matrix elements, 
denoted fy(A), play an important role, 



Im- 



1 



tn(A) 



= -v u (A),A>0. 



(2) 



Employing the relationship between the T- and S -matrices, 
Eq. (1), we can express the different v, ; in terms of phase 
shifts and the mixing angle along the physical region above 
threshold. In this way, one can write the diagonal partial waves 
as tu - (e ll5i cos2e - l)/2;p, while for the mixing amplitude 
t\2 = e' (6[+Sl) sin2e/2p. From these expressions it is straight- 
forward to obtain, for A > 0: 



v„(A)=p(A) 

y 22 (A) = p(A) 
y 12 (A) = 2p(A) 



1 - 



\ sin 2 2e 



1 



1 - cos2ecos2(5i 
4 sin 2e 



1 - cos 2e cos 2^2 
sin(<?i + 62) 
sin2e 



(3) 

(4) 
(5) 



Although not explicitly indicated, it should be understood that 
the phase shifts and mixing angle depend on A. Equation (2) 
generalizes that of an uncoupled partial wave, \mT~ x = -p, 
employed in Ref. [1]. Indeed, if we set e = in vn(A) 
and V2i(A), the uncoupled case is recovered. Note also that 
v i7 (A)/p(A) > 1. 

We apply the N/D method [2] to solve our equations for 
the T-matrix. A general NN partial wave has two types of 
cuts, the LHC and RHC, the former due to crossed channel 
dynamics and the latter to unitarity. The lightest particle that 
is exchanged between two nucleons is the pion, which deter- 
mines the onset of the LHC for A < L = -m 2 /4, with 
the pion mass. The unitarity cut occurs for A > 0. See Fig. 1, 
where the LHC and RHC are indicated separately. In the N/D 
method a partial wave fy is written as the quotient of a nu- 
merator function Njj{A) and a denominator one Djj{A). The 
function N,j only has a LHC while the function Djj has only 
a RHC. In Refs. [25, 26], a straightforward generalization of 
the one-channel N/D method of Chew and Mandelstam [2] 
was given by writing T = N ■ D 1 in matrix notation. This 
T-matrix would be symmetric, as it is required by temporal 
inversion, only under the assumption that D T (T T - T)D van- 
ishes for A — > 00 [26], where the superscript T indicates the 
transpose of the corresponding matrix. However, this is not the 
case for the chiral potentials, even at LO, e.g., in the 3 5 i- 3 £>i 
coupled partial waves. This condition is thus too restrictive 
for its application to chiral EFT, where different numbers of 
subtractions are taken in the different partial waves involved, 
whose number also varies according to the chiral order consid- 
ered in the calculation of the imaginary part of the NN partial 
wave amplitude along the LHC. 

In what follows, we generalize the procedure of Ref. [1] to 
the coupled case. Instead of making use of a matrix notation 
as in Refs. [25, 26], we write three N/D equations, one for 
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Figure 1. The thick lines correspond to the RHC and LHC, from top 
to bottom. In the same figure the integration contours Cj and Cn for 
evaluating D,j(A) and Nij(A), respectively, are shown. One has to take 
the limit e -> + . 



each of the three independent partial waves fy, as in Ref. [27], 



UA) = A 



A/A) 



(6) 



The factor A Ci > guarantees the proper threshold behavior with 
e n = I, ( 22 = €' = € + 2 and € n = {€ + C)/2 = I + 1. 
We focus here on the specific features of the coupled channel 
mechanism, referring the reader to Ref. [ 1 ] for further details 
on the general procedure followed to apply an N/D equation. 
As stated above the splitting of the fy(A) function is such that 
Nij bears the LHC and D u the RHC, and then: 

ImD, 7 (A) = -Nij(A)A^v u (A) , A > , (7) 
ImA, 7 (A) = , A < L , (8) 



with Imtij = 



Ay along the LHC. The imaginary parts of Dy 
and N/j are elsewhere along the A-real axis. 3 As argued 
in Refs. [1, 5], Ay can be calculated perturbatively in ChPT 
along the LHC, as it originates from multi-pion exchanges 
putting pion propagators on-shell. The intermediate states thus 
require at least one pion so that we apply ChPT always to irre- 
ducible A-nucleon diagrams, responsible for the discontinuity 
along the LHC. 



Two dispersion relations (DRs) can be written for the func- 
tions Djj and Ay, employing the contours C/ and Cn in Fig. 1, 
respectively. The integration along the circle at infinity van- 
ishes, if necessary, by taking sufficient number of subtractions. 
At LO in the chiral counting [5, 29], the only contribution to 
Ay along the LHC is OPE. Asymptotically, for p 2 — » -oo, 
OPE tends to constant, so that, according to the Sugawara and 
Kanazawa theorem [30, 31] one subtraction is necessary for 
the DR of Nij(A) in S wave, even though A,,(A) — > 1 /A in the 
case of OPE. On general grounds, a partial wave amplitude is 
bounded because of unitarity by constant/ y/A for A — » +oo 
so that tijDij{A)IA tij tends to constant for an S wave and 
for any other partial wave. 4 As a result the same theorem then 
requires that at least one subtraction is necessary for the S 
waves: 



D U (A) = 1 



- fdq 

x Jo 

Nij(A)=N + - fdk 

Nij(A) = - fc 
n J- 



q 2 (q 2 ~A) 

A, 7 (£ 2 )A 7 (£ 2 ) 



\dk 



k 2 (k 2 -A) 
2 \j(k 2 )D u (k 2 ) 
k 2e -J(k 2 - A) 



tu = , 



dj * . 



(9) 
(10) 
(11) 



The subtraction point is taken at threshold (see Ref. [1] for 
expressions with the subtraction point at any other position). 
One subtraction is taken for the D/j(A) function, which is fixed 
to 1 because, in view of Eq. (6), only the ratio A^y/Dy matters 
in order to determine fy. Thus, there is the freedom to fix the 
value of Djj at one point, e.g., at threshold, by simultaneously 
dividing Djj and Nij by the appropriate constant. For £y = 0, S 
wave, one subtraction is taken in A^y(A), as just discussed. In 
our present work, dedicated to the AW coupled partial waves, 
this is the case only for the 3 S i channel. The subtraction con- 
stant No is the amplitude at threshold, fn(0) = No, and then it 
can be fixed in terms of the 3 S i scattering length, a,, 



No 



Ana, 



(12) 



with the value a, = 5.424 + 0.004 fm. Below in Sec. Ill A 
we also fix No in terms of the experimental deuteron binding 
energy. 

An integral equation for the function D/j(A) results by in- 
serting Eqs. (10) or (11) into Eq. (9). However, as argued in 
detail in Ref. [1], divergent integrals appear for ( > 2 unless a 
set of €— 1 constraints is satisfied by Dy(A). These constraints 
are a generalization of those satisfied by OPE. We just quote 
the final result from Ref. [1], which is given in terms of the set 
of sum rules: 



%J —CO 



Anik^Dnik 2 ) 



k 2A 



= 0, A = 2,3,... Jij> 2 (13) 



3 Because the Schwartz reflection principle is satisfied by fy, Djj and Nu the 
discontinuity across the RHC or LHC is given by 2i times the imaginary 
part of the function. 



Here we are taking that Djj diverges as VAforA — > oo as in the uncoupled 
case [1]. This is consistent with the results obtained explicitly in this work. 
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Expanding the denominator inside the integral of Eq. (11), 
Nij(A) can be written as 



Ni,{A) 



i r L 



2 A u (k 2 )D u (k 2 ) 
k 2 (k 2 - A) 



1 r^.Anik^Duik 2 ) 




and the terms within the sum vanish if the constraints of 
Eq. (13) are fulfilled. This guarantees that N/j(A) vanishes as 
l/A CiJ , which ensures the convergence of the resulting integral 
equation for Z), ; (A). 

Let us take first £y + 0. By inserting the nonvanishing piece 
of Nij into Eq. (9), once the constraints Eq. (13) are satisfied, 
we find the following integral equation for Dy(A): 

A r L ? Aij(k 2 )Du(k 2 ) , 
Dij(A) = 1 + - dk 2 ' A ; 2 ' A J gij (A, k 2 ) , (15) 

It xj — OO 



gij (A,k 2 ) = - \dq : 
x Jo 



k 2 
Vijiq 2 ) 



(q 2 - A)(q 2 - k 2 ) 



(16) 



The functions g,j{A,k 2 ) are the generalization of g(A,k 2 ) 
given in Ref. [1] for the uncoupled case. An important techni- 
cal detail is discussed in the Appendix. We show there how the 
constraints in Eq. (13) guarantee that the functions g/j(A,k 2 ) 
are finite curing a potential divergence for ij = 22 in the 
q 2 — > limit. This divergence was noticed in Ref. [27] but 
no procedure for removing it was given there. 

The N/D method in the presence of the constraints, 
Eq. (13), was solved in Ref. [1] by means of the insertion of 
CDD poles [28] by taking advantage of the fact that the Dy 
functions are determined modulo the addition of CDD poles 
[30, 32, 33]. The main points from Ref. [1], briefly summa- 
rized, consist of using this ambiguity to include £y - 1 CDD 
poles (if tjj > 2) in Dy(A). These poles are gathered at the 
same position B, and finally the limit B — > oo is taken. The 
following equations are then obtained [ 1 ] : 



N U (A) 



n J_, 



Dij(A) = 1 + 



L ^ij{k 2 )D u (k 2 ) 
k 2e 'J(k 2 -A) 
A r^ 2 A u (k 2 )D u (k 2 ) 



(17) 



± (1 

n J-oo 



k 2 



gij(A,k 2 ) 



Z--2 

(A - B) 



e,i-\ 



(18) 



The last sum corresponds to the addition of the - 1 CDD 
poles. The coefficients c-, are determined in such a way that 
the constraints in Eq. (13) are satisfied (see Ref. [1] for further 
details). 

Note that for the P waves (€y = 1) (in the present study 
we have the mixing partial wave in the 3 S i- 3 £>i system and 
the 3 P 2 in 3 / ) 2- 3 ^ r 2 scattering), no constraints are needed [1], 
so that the sum over the CDD poles is dropped and the same 
formalism applies. This is also clear because, for this case, 
Eq. (11) vanishes as 1/A so that there is no room for restric- 
tions. 



Let us take now the case £y = 0, which only occurs for the 
3 S 1 wave. Since a subtraction is needed in Nu, Eq. (10), one 
should change Eq. (18) in two ways, as there is no sum over 
CDD poles and one has to include an extra term associated 
with the subtraction in iVy(A) for this case. It is straightfor- 
ward to obtain, by inserting Eq. (10) into Eq. (9), the appro- 
priate integral equation for Dn(A) for the 3 S \ partial wave: 



D U (A)= l-AN ogu (A,0) 



* fdV 

n J-oo 



An(k 2 )D n (k 2 ) , 2 . 
7^2 Sn(A,k ) 



(19) 



with gn(A,k 2 ) given by Eq. (16). Notice also that, from 
Eq. (10), it is clear that iVy(A) tends to constant for £y = 
and A — » oo, so that there is no need for constraints. This is 
why no sum over CDD poles is present in the previous equa- 
tion. 

To obtain the final amplitudes, the Ay(A) functions are ob- 
tained along the LHC (A < — by solving the integral 
equations in Eq. (18) or Eq. (19). Next, the functions Dy(A) 
are obtained along the RHC (A > 0) from the same equa- 
tions because the integrand is known. To obtain the functions 
Nij(A), since the constraints in Eq. (13) are obeyed, one can 
use for (/j + either Eq. ( 1 7) or the first term on the right-hand 
side of Eq. (14) (but the former is more suitable numerically, 
since it converges faster). For the 3 5 1 wave, one should use 
Eq. (10). The partial waves f,y(A) are obtained by employing 
the resulting Ay(A) and N,j(A) functions in Eq. (6). 

The main difference with respect to the uncoupled case 
treated in Ref. [1] is that now one has to solve simultane- 
ously three N/D equations for z/'=ll, 12 and 22 with the func- 
tions gij{A, k 2 ) linked between each other. They depend on the 
phase shifts 5\ and 62 and on the mixing angle e, defined in 
Eq. (1), which are also the final output of our approach. Thus, 
we employ an iterative procedure (similar to Ref. [27]) as fol- 
lows. Given an input for 6\, 62 and e, one solves the three 
integral equations for Djj(A) along the LHC, and then the am- 
plitudes for the RHC can be calculated. The new phase shifts 
<5i and 62 are obtained from the phase of the S -matrix ele- 
ments S 11 and S22, while sin2e = 2pA f|2 A^i2/|£>i2l, according 
to Eq. ( 1). In this way a new input set of functions Vy, Eqs. (3)- 
(5), results. These are used again in the integral equations, and 
the iterative procedure is finished when convergence is found 
(typically, the difference between one iteration in the three in- 
dependent D/j functions along the LHC is required to be less 
than one per mil). As initial input one can use the results given 
by UChPT [5], or some placed-by-hand phase shifts and mix- 
ing angle, and we find no dependence of our final unitary re- 
sults on the input employed. 

It can be shown straightforwardly that unitarity is fulfilled 
in our coupled channel equations, solved in the way just ex- 
plained, if |Sn(A)| 2 = |S 22 (A)| 2 = cos 2 2e for A > 0. From 
the fact that Imfi2 = vi 2 |fi2| 2 , as follows from Eq. (2), and 
sin2e = 2p\tn\ (the latter equality is valid only when conver- 
gence is reached), it follows that the phase of tu is S\ + 62, as 
required by unitarity, Eq. (1). By construction the phase shifts 
are equal to one-half the phase of the S -matrix diagonal ele- 
ments when convergence is achieved. 
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III. RESULTS 



given by 



We now present the reproduction of the phase shifts and 
mixing angles for the AW coupled partial waves with / < 3 
compared with the data from the Nijmegen partial wave anal- 
ysis (PWA) [34]. We pay special attention to the 3 5' i- 3 £>i sys- 
tem. 



E d = -^ 



(20) 



A. 3 S i - 3 D l coupled waves 

In this section we discuss our results for the 3 S \- 3 D\ cou- 
pled waves. Previous papers applying the N/D method to ad- 
just AW scattering are Refs. [21-24]. We already commented 
about Ref. [21, 24]. The other two works, by Wong and Scotti 
[22, 23], include, together with OPE, other heavier mesons, 77, 
p, to, and <p is also included in Ref. [23]. Thus, these works 
follow the basic ideas of meson theory of nuclear forces, 
that were also used for the construction of AW potentials 
[35]. There are some approximations in Refs. [21-23] that 
we avoid in our work. For example, only elastic unitarity is 
used in Refs. [21, 22] neglecting the mixing between coupled 
partial waves. Reference [23] considers the mixing only for 
3 S\— 3 D\ coupled partial waves. In addition, in order to sat- 
isfy the threshold behavior for partial waves with I > 2, so 
that they vanish as A 1 , Refs. [22, 23] make use of a rather 
ad-hoc formula. This method was criticized in Ref. [36] be- 
cause it includes an unphysical pole for every partial wave at 
a CM squared energy si, somewhat below Am 2 (the thresh- 
old for AW scattering). In addition, Refs. [22, 23] also have a 
cutoff dependence in the way the vector resonance exchanges 
damp to avoid their divergences at infinity. Though the results 
of Refs. [21-23] are interesting and typically obtain a good 
reproduction of data at the phenomenological level, we offer 
here a novel way of employing the N/D method in light of 
EFT. We then present the method ready to be used in a sys- 
tematic way by improving, order by order, the discontinuity of 
the partial wave amplitudes along the LHC as it involves only 
AW irreducible diagrams, as discussed above [5]. We satisfy 
exact unitarity for all the partial waves as well. It is also im- 
portant to stress that the N/D method for coupled channels 
is now presented in a way ready to be used at any chiral or- 
der, without being constrained to satisfy the too demanding 
Bjorken-Nauenberg condition [26] in order to end with sym- 
metric partial waves. We accomplish the right threshold be- 
havior for t > 2 by adding CDD poles at infinity, which is 
always legitimate in the N/D method if there are good rea- 
sons to include them (which have been offered before [1]). 
Thus, we do not need to modify the right analytical properties 
of partial waves by including a fictitious pole in s\ which is 
then fine tuned to data, as done in Refs. [22, 23]. 

The deuteron (d) is a neutron-proton (np) bound state with 
total angular momentum / = 1 and spin 5=1 (and isospin 
0). As such, it is seen as a pole below threshold (|p| 2 < 0) in 
the physical Riemann sheet in 3 S i- 3 £>i coupled partial waves. 
The binding energy of the deuteron, Ed (defined positive), is 



where k 2 d is the three-momentum squared at which the pole 
is located, so that it is negative. Specifically, in our approach 
it appears as a zero in the functions D,j(A). From the ampli- 
tudes calculated in Sec. II we find the deuteron at the posi- 
tion k d = -0.08m 2 m tne amplitude, corresponding to 
Ed 1.7 MeV. Recall that the subtraction constant A^o ap- 
pearing in the 3 S 1 partial wave is determined by fixing the 
3 S 1 scattering length to its experimental value, Eq. (12). There 
is still a remnant input dependence for the 3 S \- 3 Di coupled 
partial waves in our unitary solutions that we fix by requir- 
ing that the deuteron pole position is the same in the 3 S 1 and 
in the mixing partial wave. Independently of the input we do 
not find any pole in the 3 D\ partial wave. Indeed, if we dis- 
regard the coupling between 3 5 1 and 3 D\ and use the method 
of Ref. [1] for uncoupled waves, the pole appears in the same 
position in 3 5i and, again, it does not appear in 3 £>i. How- 
ever, the pole should be located at the same energy in every 
channel, but this is not the case because we are not using a 
matrix formalism but solving the three linked N / D equations 
independently. Notice also that the deuteron is found mainly 
in a 3 5 1 state, and thus the coupling to 3 D\ is very weak. 

In order to cure this deficiency and having the right pole 
structure guaranteeing the presence of the deuteron pole in 
3 Z?i, we write a twice subtracted DR for the 3 D\ partial wave, 
such that the function DniA) has a zero at a given k 2 d . The DR 
reads 



Du(A) = 1 



*3 



A(A-k 2 d ) f+~ Vijiq^NtM^q 21 " 



^5) r 

x Jo 



q\q 2 - A)(q 2 - k 2 ) 



written in a way that is valid both for the 3 D\ partial wave 
(ij = 22) and for the mixing partial wave (ij = 12), although 
we do not use it for the latter. By inserting the expression for 
Nij(A), Eq. (11), into the previous equation, we end up with 
the following integral equation 



A/A) = 1 - - 2 



A(A - kj) T L ;2 A u (k 2 )D u (k 2 ) 



^ fdV 

n J-co 



g { "\A,k 2 ), (22) 



k 2e,j *u 



where g\f(A, k 2 ) is a generalization of the functions gij(A, k") 
inEq. (16), 



(J) 2 1 f + % v n {q 2 )q 2{t '^ 



71 Jo ' (q 2 -A)(q 2 -k 2 )(q 2 -k 2 d ) 



(23) 



For 3 S \- 3 Di waves, we have I — and £' = 2, so that £12 = 1 
and ^22 = 2, and the previous integrals are convergent because 
of the extra subtraction taken. Recall that, in the formalism 
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Figure 2. (Color online) Comparison of our results for the 3 S 1 and 3 Di phase shifts and the mixing angle e\ to the Nijmegen PWA [34] , shown 
by the dot-dashed (red) lines. The solid (black) lines correspond to fixing the 3 S i scattering length to experiment, while the dashed (blue) lines 
in addition fix the deuteron pole position in 3 D\ at the same value as that in 3 5i. On the other hand, the double dotted (green) lines stem by 
fixing the deuteron pole position in the 3 S \ partial wave at its experimental value. The dash-double-dotted (cyan) lines correspond to having 
additionally fixed the deuteron pole in the 3 D\ partial wave at the same point as in 3 Si. 



first presented in Sec. II, one must take into account a con- 
straint for the D22(A) partial wave in order to end with a con- 
vergent integral equation. Note that from Eqs. (22) and (23) 
the high-energy behavior of the functions Dy changes, now di- 
verging as A 3 / 2 , instead of A 112 as in Sec. II or in Ref. [1]. As a 
result, the criterion of imposing that Nij — > 1 /A '•' for A — » oo, 
the one used in Ref. [1] to deduce the need of constraints, 
does not hold in this case because of the extra subtraction. 5 
The price to pay for having included the second subtraction is 
the need for an input value for k 2 p which has to be provided. 
It is then more natural for the 3 S i- 3 £>i system to fix the bind- 
ing energy of the deuteron to its experimental value than the 
scattering length, as we do below. 

As stated in Sec. II, an iterative procedure is followed in 
order to obtain our final results for the phase shifts and the 
mixing angle from the three N/D equations coupled. For ev- 
ery iteration along that procedure, one obtains from the 3 5' i 
wave amplitude the deuteron pole position, ki. This is the 
value used as an input for the function £>22(A) at every step. 
In this way it is not fitted as a free parameter in order to fix 
the deuteron binding energy, but it comes out in a natural 
way from 3 5 1 and the coupled-channel mechanism. The re- 
sults that we obtain with this approach are shown in Fig. 2 
by the dashed (blue) lines, while those obtained when there 
is no deuteron pole in 3 Di, using Eq. (18) instead of Eq. (22) 
with Cjj = 2, correspond to the solid (black) lines. The re- 
sults are compared with the Nijmegen PWA [34] given by 
the dash-dotted (red) lines. For the 3 5' i phase shifts both lines 
are very similar. The differences are larger for the 3 £>i phase 
shifts, which are then quite sensitive to reproducing correctly 
the deuteron pole also in the 3 £>i partial wave. Indeed, the re- 
sult without imposing the deuteron in this partial wave is very 
similar to that obtained from perturbative OPE [37]. Differ- 
ences are rather small for the mixing angle e\ . As the main 
contribution to the deuteron comes from 3 5 1, its position re- 



5 From Eq. (11) it follows immediately that N22 — > 1/A which is the be- 
havior required for N22 - t2iD 22 1 A 2 , taking into account the high-energy 
behavior of 022(A) just discussed. 



mains almost unchanged compared with the uncoupled case, 
with a value obtained for the binding Ed 1 .7 MeV, once the 
experimental scattering length is fixed. This corresponds to an 
effective range r 0.46 fm, which is much smaller than the 
experimental value r — 1 .749 fm, the difference being around 
a 70%. This fact is already well documented in the literature 
[33]. Indeed, Ref. [21] shows that when the N/D method is 
used with only OPE as the source of the imaginary part along 
the LHC, one needs to fit two experimental inputs for every 
AW S wave in order to reproduce the scattering length and ef- 
fective range. For 3 S 1 the scattering length and the deuteron 
binding energy are taken (we take the same input in Sec. Ill B 
below), while for 'So two well measured phase shifts at differ- 
ent energies are employed. This result from Ref. [21], and our 
own ones presented below in Sec. IIIB, makes us confidence 
that a NLO study in ChPT with the N/D method will be phe- 
nomenologically successful because a new counterterm enters 
at this order, multiplying an energy dependent monomial. The 
authors of Ref. [21] make the approximation of considering 
only elastic unitarity for 3 5 1 , neglecting its coupling with 3 Dj , 
while our treatment is exact. 

It is also interesting to fix the subtraction constant iVo in 
terms of the deuteron binding energy and then compare with 
our previous results when the scattering length was fixed. Im- 
posing D\\{k 2 d ) = from Eq. (19) and solving for No, one has 

1 + § fdk 2 A u (k 2 )d u (k 2 )gu(k 2 ,k 2 )/k 2 

No = -r , (24) 

k 2 ( 8u (k 2 p 0) + g(k 2 d j) 

where we have first split 

D„(A) = £/ 11 (A)-* 2 iVbgii(A,0), (25) 

from where the function d\\{A) is defined. We have also intro- 
duced in Eq. (24) the function @(A) given by 

1 C L 

§{A) = - \dk 2 A n (k 2 )gn(k 2 ,0)gn(A,k 2 ). (26) 
n J-oo 

The integral equation for d\\{A) can be obtained from that in 
Eq. (19) taking into account Eq. (25) and replacing No with its 
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expression Eq. (24). This results in 

dn(A) = 1 + - Jdk ^ ig u (A,k) (27) 

g n (k 2 ,k 2 d )g(A) \ A @{A) 
gu(k 2 d ,0) + g(fyj k 2 d gu(k 2 d ,0) + g(k% ' 

As in the previous case we fix the dependence on the in- 
put by requiring that the deuteron pole in the mixing wave is 
located at the same position as in the 3 S i wave, at the k 2 d corre- 
sponding to the binding energy Ed - 2.2 MeV. Regarding the 
3 Di partial wave no pole position is found unless one imposes 
it in the £>22(A) function, making use of Eq. (22), having then 
the right pole structure. Once the deuteron pole is imposed the 
value that we obtain for the 3 S i scattering length is 4.6 fm and 
for the effective range 0.41 fm. The latter is indeed very sim- 
ilar to the values obtained before when the scattering length 
was taken as input. The resulting scattering length is about 
15% lower than its experimental value. We show in Fig. 2 by 
the dash-double-dotted (cyan) lines the results obtained when 
the deuteron pole is imposed in the 3 S i and 3 D\ partial waves, 
while the double-dotted (green) line is for the case when the 
deuteron pole position is imposed only in the former. The re- 
sults are rather similar to the case when the scattering length 
was fixed. The most sensitive observable is the mixing angle 
e\ where the largest difference happens in the peak, somewhat 
less than 1°. 

It is worth comparing our results with the pionless effective 
field theory. In this case pions are integrated out as heavy de- 
grees of freedom. We can reach this limit by taking gA — > in 
our results, which implies A, ; = 0. Only the term proportional 
to No survives in Eq. (19) and N n (A) = No from Eq. (10). 
We can determine No by fixing the experimental scattering 
length, Eq. (12), or by reproducing the deuteron binding en- 
ergy No = —An I s/m 3 E l i- The former case is given by the solid 
(black) line and the latter by the double-dotted (green) one in 
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Figure 3. (Color online) Comparison between the results obtained for 
the theory without pions and our full results at LO. The Nijmegen 
PWA [34] corresponds to the dot-dashed (red) lines. When the 3 5 1 
scattering length is fixed one has the dashed (blue) line for the pio- 
nfull case and the solid (black) line for the pionless one. When the 
deuteron binding energy is fixed the dash-double-dotted (cyan) line 
for the pionfull theory and the double-dotted (green) line for the pio- 
nless case result. 



Fig. 3. For comparison we also show the lines corresponding 
to our full results, obtained by fixing the scattering length and 
the deuteron binding energy to their experimental values. The 
former case corresponds to the dashed (blue) line and the latter 
to the dash-double-dotted (cyan) line, as shown in Fig. 2. One 
observes that the inclusion of pions significantly improves the 
phase shifts and also makes the results more stable indepen- 
dently of whether the scattering length or the deuteron pole 
are adjusted. 



B. One extra subtraction 

Now we impose that the 3 5' i partial wave reproduces the 
experimental values for the 3 S i scattering length, a,, and the 
deuteron binding energy simultaneously. Similar restrictions 
were already considered in Refs. [21, 24]. To accomplish it we 
introduce one extra subtraction constant in the Dn(A) func- 
tion by taking one more subtraction in the DR. In this way we 
enhance the role played by the low-energy region because the 
extra subtraction gives more weight to the low-energy part of 
the integrand in the DR, so that it vanishes more rapidly as 
A — > oo. The new DRs for Nn(A) and D n (A) read 



Nu(A) = N + 
D n (A) = l 



A f L d! _ 2 An(k 2 )Dn(k 2 ) 

n J_ M k 2 (k 2 -A) ' 

A A(A-k 2 .) f°° 

d -N dk 2 

n Jo 



vi i(k 2 ) 



k 2 



(k 2 - A)(k 2 - k 2 )k 2 
A(A-k 2 d ) r L ji2 A u (k 2 )D u (k 2 ) 



n J_ m k L 



with 



gif (A, k 2 ) = - fV 
n Jo 



vn(q 2 ) 



(q 2 - k 2 )(q 2 - A)(q 2 - k 2 ) 



(28) 



(29) 



By construction D\\(k 2 d ) — in Eq. (28), which guarantees the 
presence of the deuteron in its experimental position. Having 
the right value for the 3 S \ scattering length fixes the constant 
A^o to Eq. (12). The extra subtraction taken in D\\, Eq. (28), 
will also be studied when considering the NLO ChPT contri- 
bution to the discontinuity across the LHC because then the 
resulting A,y(A) diverges as A for A — > oo. 

The deuteron pole is also imposed in the 3 D\ partial wave 
by employing Eq. (22) so that the right pole structure is ac- 
complished. The input is fixed such that the resulting deuteron 
pole position in the mixing partial wave 3 5' i- 3 £>i is located in 
the same position as for the other two coupled partial waves, 
as already discussed above. 

In Fig. 4 we show, from left to right, the 3 S i and 3 D\ phase 
shifts and the mixing angle e\ resulting from Eq. (28), in that 
order. A clear improvement as compared with Fig. 2 is ob- 
served, so that now the resulting curve run closer to the Ni- 
jmegen PWA [34] for the 3 S i phase shifts. An improvement 
also happens for the mixing angle e\ which now overlaps bet- 
ter with the Nijmegen results for three-momentum up to about 
100 MeV and later the trend of the curve tends to follow that 
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Figure 4. (Color online) From left to right we plot our results, solid (black) curves, for the 3 5i and 3 Di phase shifts and the mixing angle 
e,, when the experimental 3 S \ scattering length and deuteron binding energy are imposed. The Nijmegen PWA [34] data are shown by the 
dot-dashed (red) lines. 



of the Nijmegen PWA. Let us also stress that the failure to 
reproduce e\ in the Kaplan-Savage-Wise scheme [14] was the 
main reason to conclude that its perturbative treatment of pion 
exchange was not appropriate [38]. In contrast, our LO re- 
production of e\ in Fig. (2) is already quite close to the Ni- 
jmegen results [34] and improves when considering the extra 
subtraction, as shown in Fig. 4. This is a clear indication that 
e\ will be also properly reproduced at NLO in the calcula- 
tion of Ay(A), although the adjusted value for this subtraction 
constant would change due to the addition of the two-pion ex- 
change contributions. 

Next we evaluate the three independent deuteron parame- 
ters that can be calculated from AW scattering [39]. The first 
quantity is the binding energy of the deuteron that is fixed to 
its experimental value as input. The second quantity that we 
consider is the asymptotic D/S ratio 77. For that we make use 
of the Blatt and Beidenharn parameterization [40] and diago- 
nalize the 3 S \- 3 Di S -matrix, S i, by an orthogonal real matrix 
O, 

Si = OS l;diag<5 1 , 




cos e\ — sin e\ 
sin e\ cos e\ 



S^^[ S ° S2 ). OO) 

In terms of e\ one can write for the asymptotic D/S ratio 77 as 
[39,41] 

77 = - tan €\ . (31) 
The third quantity that we calculate is ; times the residue of 



the eigenvalue So at the deuteron pole position a = 




So — ^— : + regular terms . (32) 

a + ;|p| 

We should remark that because we do not employ AW po- 
tential to study AW scattering we cannot compute the wave 
function of the deuteron and in terms of it evaluate straight- 
forwardly (in the simplest approximation) other quantities, 



e.g., the deuteron electric quadrupole moment Q or the mean- 
square deuteron radius (r 2 ) 1 ^ 2 . This does not mean that we 
cannot obtain such observable quantities from our T-matrix 
but simply that we should consider other processes beyond 
pure AW scattering. For instance, in order to calculate the 
mean-square deuteron radius (r 2 ) we should proceed as in 
Ref. [42] to calculate the same quantity but for the /o(500) or 
cr resonance, where nn scattering in the presence of a scalar 
source was calculated. Similarly, we should study here AW 
scattering in the presence of a scalar source giving rise to the 
matter form factor of the deuteron. This is beyond the present 
study and requires an independent study. 
The resulting values that we obtain are 

77 = 0.028 , N 2 p = 0.74 frrT 1 . (33) 

Our results compare well with the experimental determina- 
tions, 77 = 0.0271(4) [43] and 77 = 0.0263(13) [44]. They are 
also close to those evaluated in Nijmegen PWA 1993 [34] 

77 = 0.02543(7) , N 2 p = 0.7830(7) frrT 1 . (34) 

Thus, once we reproduce simultaneously the deuteron bind- 
ing energy and the 3 S 1 scattering length, the other properties 
of the deuteron that can be extracted from scattering compare 
well with the values determined in partial-wave analyses or 
experiment. 

We obtain the following value for the effective range r, 

r= 1.56(3) fm, (35) 

where the error is just statistical by fitting the low-energy 
phase shifts generated by our own amplitudes. This num- 
ber is quite close to the Nijmegen PWA 1993 [34] result, 
r = 1.753(2) fm. 

In Ref. [45] the OPE potential from ChPT is employed 
in a LS equation solved by making use of an interesting 
method based on identifying the input with the T-matrix deep 
in the LHC, writing the potential in terms of it. Their re- 
sults for 77 and r are very similar to ours in Eqs. (34) and 
(35), obtaining the intervals of values 77 = 0.0281-0.0293 
and r = 1.36-1.58 fm. Their results for the elastic 3 S 1 phase 
shifts are also quite similar to ours, though for 3 D\ they are 
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closer to Nijmegen points [34]. Regarding the mixing angle 
ei, Ref. [45] obtains that for a large renormalization scale // 
the resulting curves depart from the Nijmegen data [34] by 
an absolute amount similar to ours for |p| > 100 MeV (our 
results lie above, while theirs lie below). One should keep in 
mind that we have taken the scattering length and the bind- 
ing energy as input for our calculations, while Ref. [45] only 
adjusts the scattering length. 

It is well known since the 1960s that for 3 S i- 3 £>i coupled 
partial waves, solving a LS equation in terms of the OPE po- 
tential gives a significantly better phenomenology than solv- 
ing the N/D method taking for Ay(A) the discontinuity along 
the LHC induced by OPE [36]. However, it is worth keep- 
ing in mind that [45], as well as [10], obtain phase shifts for 
'So which are very similar to ours in Ref. [1]. It is known 
that the 'So phase shift data of Nijmegen [34] are reproduced 
quite closely [9] once two-pion exchange contributions and 
NLO LECs in the four-nucleon Lagrangian are included. In 
our novel theory, which calculates the AW partial waves from 
ChPT by employing the N/D method, there is no reason to 
expect that the phase shifts should be reproduced at LO worse 
in the 'So partial wave than in the 3 Si- 3 Di coupled waves 
(which results strictly correspond to Fig. 2 in terms of only 
one subtraction being needed). In this respect, it is reward- 
ing that, by considering NLO contributions to the AW poten- 
tial in the standard Weinberg approach [9], one can obtain 
good results for 'So- This should be also expected for the 
3 Si- 3 Di case within our approach. Indeed, we have already 
seen that, by including one extra subtraction, the reproduction 
of phase shifts (particularly for 3 S i) and mixing angle clearly 
improves. When considering two-pion exchange at NLO some 
extra counterterms are needed because A,j(A) diverges as A for 
A — » -oo along the LHC. 

Solving a LS equation with OPE for the 3 S i- 3 £>i system 
is much more successful phenomenologically than for the 
'So case. One should be aware that this is something that is 
checked a posteriori and is not rooted in the chiral counting 
(in which our approach is based). From our point of view the 
ladder resummation in the LS for the 3 S i- 3 £>i case provides 
higher orders terms to A,y(A) in the right direction. However, 
this improvement should come out when applying the N/D 
method to (just a few) higher orders, because along the LHC 
Ajj(A) is perturbative and amenable to a chiral expansion as 
discussed. For 'So the higher orders in A(A) provided by the 
LS equation are not the important source of dynamics and one 
has to really consider the full machinery in order to incorpo- 
rate at higher orders two-pion exchange with the associated 
chiral counterterms. It is our aim to develop for the time being 
a NLO (or NNLO) study of AW scattering with our approach 
based on the N/D method and the ChPT calculation of A(A) 
in order to definitively settle this important issue. We would 
like to stress that at this stage our study is mostly exploratory 
and not competitive with the current sophisticated potentials 
[34] or calculated at higher orders from ChPT [8, 9]. 

The set of works [19, 46-49] gives rise to a remarkable de- 
scription of deuteron properties employing the AW potential 
given by OPE in a LS equation, e.g., Ref. [19] achieves for 
many observables a 2 - 3% of deviation with respect to the 
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Figure 5. (Color online) Comparison of our results for the 3 P 2 , 3 F 2 , 
3 Z) 3 and 3 Gt, phase shifts and the mixing angles e 2 and e 3 , shown by 
the solid (black) lines, to the Nijmegen PWA [34], represented by the 
dot-dashed (red) lines. 



experimental values. But this is not the only aim of an EFT. 
That is, one does not expect such a high degree of convergence 
by taking only the LO ChPT AW potential. This is more a mat- 
ter of phenomenological success and not rooted in the chiral 
EFT. For baryon ChPT the expansion scale is not so great, 
A a l27T 2 fx/g\tn a 500 MeV [5, 50], and such a great preci- 
sion is thus difficult to understand from the ChPT expansion. 
We want to emphasize this point (consider, e.g., the not so 
great achievement for the 'So case) and develop a formalism 
where contributions to a given process can be obtained order 
by order systematically in the chiral EFT expansion of A(A). 



C. Higher partial waves 

In this section, we present the results for the spin triplet 
waves with total angular momentum J — 2 and 3, obtained 
with the formalism derived in Sec. II. They are shown by the 
solid (black) lines in Fig. 5, where they are compared with the 
Nijmegen PWA [34] [dash-dotted (red) lines]. 

We already see a good agreement with data for ^Fj and 3 G3 
as well as for the mixing angles €2 and 63. The lower partial 
waves 3 P2 and 3 £>3 are not well reproduced with only OPE 
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yet. This fact for the 3 £>3 partial wave was already observed 
in Ref. [37], where OPE was treated perturbatively. In this 
reference 3 D 3 is also obtained with opposite sign to the data. 
In Ref. [10], with one counterterm promoted to LO for the 3 P2 
wave, the situation is similar. The 3 P2 and 3 £>3 phase shifts are 
not well reproduced at LO, while the others compare well with 
the data. We expect to restore the agreement with experiment 
at higher orders in the application of our method to 3 Pi and 
'D 3 . 



IV. CONCLUSIONS 

We have developed a new set of equations for the N/D 
method in coupled partial waves, extending our previous 
work, Ref. [1], restricted to uncoupled partial waves. This 
method is presented in a novel way, adequate to improve the 
results systematically by taking higher orders in the chiral ex- 
pansion of the calculation of the discontinuity of the partial 
wave amplitudes along the LHC, A, ; (A) with A < 4m 2 . This 
extension is accomplished by providing three N/D equations 
for each set of partial waves coupled. The solution is obtained 
in an iterative and self-consistent way. The correct solution 
satisfies unitarity and for the case of the 3 Si- 3 Di system the 
deuteron pole is located at the same position in all the waves, 
having the correct pole structure. 

As in Ref. [1] our approach guarantees the right thresh- 
old behavior for partial waves with orbital angular momen- 
tum iij > 2 by satisfying £y — 1 constraints. Since the function 
Djj(A) is determined modulo the addition of Castillejo-Dalitz- 
Dyson poles (that correspond to zeros of the AW partial waves 
along the real axis) we have then added £y - 1 of such poles 
at infinity in D,j for £y ^ 2. By sending such poles to infinity 
we do not include any zero of any AW partial wave at finite 
energies. In addition, the residues of these poles in D,j{A) are 
fixed once the sum rules are satisfied, so that no new parame- 
ters are included. At low energies the CDD poles behave like 
adding a polynomial of degree tjj - 2 to Dy(A). 

We have studied the 3 Si- 3 Di coupled waves by fixing the 
resulting subtraction constant to either the experimental value 
of the 3 S 1 scattering length or the deuteron binding energy. 
We find that the 3 D[ phase shifts are the most sensitive to this 
choice. As expected, in all cases the triplet S wave effective 
range comes out much smaller than in experiments. Then, we 
added one extra subtraction to calculate the 3 S i wave requir- 
ing the simultaneous reproduction of the deuteron binding en- 
ergy and triplet S wave scattering length. The resulting 3 S 1 
phase shifts are much improved and we then obtain the effec- 
tive range and deuteron properties close to their experimental 
values. We have also considered the pionless case and com- 
pared with our full results, which include OPE. It was seen 
then that the results clearly improve for the latter case. For 
the waves with orbital angular momentum t 5= 1 at LO there 
is no subtraction constant and the results are parameter free. 
The resulting phase shifts and mixing angles agree well with 
the Nijmegen PWA results, except for the 3 P2 and 3 £>3 partial 
waves. 

Certainly, including OPE as the only source of discontinu- 



ity along the LHC is phenomenologically just a first step and 
a NLO calculation should be undergone to establish the capa- 
bility of the method to reproduce properly AW scattering data. 
However, one should stress at this point that our approach 
based on the N/D method offers a way to calculate AW scat- 
tering independently of cutoff, because only convergent inte- 
grals appear, while keeping the chiral power counting. The 
dispersive integrals are convergent by taking the appropriate 
number of subtractions with the related subtraction constants 
fixed to experimental data. At LO only two subtraction con- 
stants appear in the 'So and 3 S i partial waves, the same num- 
ber as LO ChPT counterterms [6]. This method allows one to 
perform calculations systematically, order by order, in ChPT. 
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Appendix: The gij(A, k 1 ) functions 

On general grounds, the following threshold behavior are 
found for the v,j(A) functions: 



vn(A)ocA 1/2 vi 2 (A)ocA 



-1/2 



v 22 (A) oc A 



-3/2 



This can also be seen by inserting the low energy behavior 
of 6i (oc A M/2 ), 6 2 (oc A r+1/2 ), and e (oc A {M ' +1)/2 ) in the ex- 
plicit expression for v,y(A), Eqs. (3)-(5). No problem occurs in 
the integrand for the functions g\i(A, k 2 ) and gn(A, k 2 ), when 
these low-energy behaviors are inserted, but the divergence 
in V22(A) could lead to a divergence in the function g2i(A, k 2 ). 
This was already pointed out in Ref. [27], as a potential source 
of divergences. However, a more careful analysis shows that 
this divergence vanishes owing to the sum rules, Eqs. ( 1 3). For 
the g22(A, k 2 ) integral one has 



S22(A,£ 2 ) 



1 r + °° 

* Ja->o 



2v 



nAk 2 VI 



(q 2 - A)(q 2 - k 2 ) 
+ regular terms 



where v 2 2(A) = voA -3 ^ 2 for A — > 0. In the previous equation 
the regular terms refer to the rest of the contributions to the in- 
tegral, which do not diverge for A — > 0. The divergent term in 
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the previous equation enters into the integral Eq. (18) through 
the function g2i(A, k 2 ), giving rise to a term proportional to 



^ 2 A 22 (ft 2 )P 22 (fc 2 ) = 



which vanishes owing to the constraints of Eq. (13). Note that 
every channel in which g2i(A) is involved has I' > 2 (the 
lowest value for (' corresponds to the 3 D\ wave), and thus the 
sum rule above applies. The constraints Eq. (13) thus show a 
new important facet beyond the original motivation for their 
introduction. 
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